Precision Rosenbluth measurement of the proton elastic form factors 
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We report the results of a new Rosenbluth measurement of the proton form factors at Q val- 
ues of 2.64, 3.20 and 4.10 GeV 2 . Cross sections were determined by detecting the recoiling proton 
in contrast to previous measurements in which the scattered electron was detected. At each Q 2 , 
relative cross sections were determined to better than 1%. The measurement focussed on the ex- 
traction of Ge/Gm which was determined to 4-8% and found to approximate form factor scaling, 
i.e. {i p Ge ~ Gm. These results are consistent with and much more precise than previous Rosen- 
bluth extractions. However, they are inconsistent with recent polarization transfer measurements 
of comparable precision, implying a systematic difference between the two techniques. 

PACS numbers: 25.30.Bf, 13.40.Gp, 14.20.Dh 



Reproducing the structure of the proton is one of the 
defining problems of QCD. The electromagnetic structure 
can be expressed in terms of the electric and magnetic 
form factors, Ge and Gm, which depend only on the 
4-momentum transfer squared, Q 2 . They have tradition- 
ally been determined utilizing the Rosenbluth formula Q 
for elastic e-p scattering: 
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where r = Q 2 /4M 2 , e is the longitudinal polariza- 
tion of the exchanged virtual photon, = 1 + 2(1 + 
r)tan 2 (8 e /2), M p is the proton mass, and e is the elec- 
tron scattering angle. The form factors are related to the 
spatial distributions of the charge (Ge) and magnetiza- 
tion (Gm) in the proton, and in the non-relativistic limit 
are simply the Fourier transformations of these distribu- 
tions. 



A Rosenbluth separation is performed by varying the 
incident electron energy and electron scattering angle to 
keep Q 2 constant while varying e. Rosenbluth separa- 
tions of Ge and Gm have been reported from 1960 to 
the present day (see Refs. 0,0 and references therein). 
Fits to these data yield ^ p Ge/Gm ~ 1 implying 
similar charge and magnetization distributions. At large 
Q 2 values, Gm dominates the cross section at all e val- 
ues (contributing more than 90% for Q 2 > 4 GeV 2 ), and 
thus while a Rosenbluth separation can yield a precise 
extraction of Gm, the uncertainty in Ge increases with 
increasing Q 2 . 

The high-Q 2 behavior of the electric form factor can be 
more precisely determined in polarization transfer exper- 
iments, where longitudinally polarized electrons are scat- 
tered from unpolarized protons and both transverse and 
longitudinal polarization are transferred to the struck 
proton. During the last few years, polarization transfer 
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experiments have been performed at Jefferson Lab which 
measured G E /G M up to Q 2 = 5.6 GeV 2 HQ. These 
measurements show the ratio decreasing with increasing 
Q 2 in stark contrast to the approximate scaling observed 
in Rosenbluth measurements. 

At high Q 2 , the quoted uncertainties on the polariza- 
tion transfer results are much smaller than those for the 
Rosenbluth extractions. This fact, combined with the 
scatter in the results from different Rosenbluth measure- 
ments, led to speculation that the Rosenbluth determi- 
nations oIGe were unreliable. While the scatter appears 
to be the result of incomplete treatment of normalization 
uncertainties when combining data from different exper- 
iments Q, the Rosenbluth technique is still very sensi- 
tive to small corrections to the cross section at large Q 2 
values. In this letter, we report the results of a new ex- 
periment that utilizes an improved Rosenbluth technique 
to determine the proton form factor ratio Ge/Gm with 
uncertainties a factor of two to three smaller than all pre- 
vious Rosenbluth measurements and comparable to the 
precision of the polarization transfer measurements. 

Experiment E01-001 was performed in Hall A at Jeffer- 
son Lab. A 70 /xA electron beam with energies from 1.9 to 
4.7 GeV impinged on a 4-cm liquid hydrogen (LH2) tar- 
get. Protons from elastic e—p scattering were detected in 
the High Resolution Spectrometer. An aerogel Cerenkov 
detector was used to eliminate charged pions from the 
data, and data from a "dummy" target were used to 
measure the contribution from the aluminum walls of the 
target. Table [I] lists the kinematics of the experiment. 
Detailed descriptions of the spectrometer and beamline 
instrumentation can be found in Ref. [7J. 



TABLE I: Kinematics of the experiment. 
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Because the fractional contribution of Ge to the cross 
section is small at large Q 2 , its extraction is highly sen- 
sitive to corrections that modify the e dependence. All 
previous Rosenbluth separations involved detection of the 
electron. The rapid variation of the electron cross sec- 
tion with scattering angle means that any rate-dependent 
corrections to the cross section will have a strong e de- 
pendence. Varying e also leads to large changes in the 
momentum of the scattered electron, and so momentum- 
dependent corrections will affect the extracted value of 
Ge- Finally, radiative corrections to the cross section 
also yield significant e-dependent corrections. 
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FIG. 1: (Color online) The measured Sp spectrum for the low 
(top) and high (bottom) e points at Q 2 = 3.2 GeV 2 (circles). 
The dotted (magenta) line is the background from the target 
walls, the long-dash line (green) is the simulated background 
from 7p — + ir°p and jp — > jp reactions, the short-dash (blue) 
line is the simulated elastic spectrum, and the solid (red) line 
that goes through the data is the sum of the target wall, 
elastic, and background contributions after each contribution 
is normalized to the data (see text). 



Detecting the struck proton greatly reduces all of these 
sources of uncertainty. The proton differential cross sec- 
tion (da /dfl p ) changes by less than a factor of two over 
the e range of the measurement while the electron cross 
section (da/dtt e ) varies by almost two orders of magni- 
tude. In addition, the minimum cross section is twenty 
times larger for the proton than for the electron. The 
proton momentum is constant at fixed Q 2 , while the cor- 
responding electron momentum varies by a factor of ten. 
Additionally, the corrections due to beam energy offsets, 
scattering angle offsets, and radiative corrections (in par- 
ticular the electron bremsstrahlung corrections) all have 
a smaller e dependence when the proton is detected. 

Figure ^ shows the "missing momentum" spectra at 
two e values for Q 2 = 3.20 GeV 2 . The missing momen- 
tum, Sp, is defined as the difference between the mea- 
sured proton momentum and the proton momentum ex- 
pected for elastic scattering at the measured proton an- 
gle. The Sp spectrum is dominated by the e-p elastic 
peak at Sp » 0. To precisely separate elastic scatter- 
ing from other processes, we need to know the shape of 
the elastic peak. The peak is simulated using the Monte 
Carlo code SIMC which takes into account the accep- 
tance and resolution of the spectrometer, energy loss and 
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multiple scattering of the proton, and radiative correc- 
tions |8| • The resolution of the simulation has been mod- 
ified to reproduce the small non-gaussian tails observed 
in the data. These are matched to the coincidence data, 
taken for two beam energies at Q 2 —2.64 GeV 2 , where the 
background contributions that dominate the singles spec- 
trum at large \Sp\ values are strongly suppressed. The 
resolution of the elastic peak is dominated by the angu- 
lar resolution. The greater width at large e is due to the 
increased sensitivity of proton momentum to scattering 
angle. 

Also shown in Figgis the decomposition of the back- 
ground into two components. The background that ex- 
tends to high Sp is due to quasielastic scattering and 
other reactions in the target walls. The other background 
is mainly due to jp — > tt°p events, with a small (1-4%) 
contribution from jp — > jp. The spectrum from these re- 
actions was modeled using a calculated bremsstrahlung 
spectrum and an s~ 7 cross section dependence. Because 
of the finite pion mass, the proton spectrum from pion 
photoproduction cuts off approximately 10 MeV below 
the elastic peak. For small proton angles, where we have 
the best resolution, the background from pion photopro- 
duction can be cleanly separated from the elastic protons 
while at large angles the background becomes small. The 
correction to the elastic cross section due to contributions 
from the target walls is approximately 10%, while the in- 
elastic processes from hydrogen contribute less than 2%. 

Because the thicknesses are different for the LH2 
and dummy targets, the bremsstrahlung yields are also 
slightly different. We use the dummy data to determine 
the shape of the endcap contributions, but normalize the 
contribution to the LH2 spectrum at large Sp, where the 
contribution from the hydrogen is negligible. While the 
shape of the bremsstrahlung spectrum differs slightly be- 
tween the dummy and LH2 targets, the effect is only 
noticeable near the endpoint. The small uncertainty due 
to the difference in shape is taken into account in the 
systematic uncertainties. 

After removing the background due to the endcaps, 
the simulated ir° photoproduction spectra were normal- 
ized to the low-momentum sides of the Sp spectra (taking 
into account the elastic radiative tail). Removing this 
background yields clean spectra of elastic events which 
are then compared to the elastic simulation. The elastic 
cross section is taken to be the value used in the simula- 
tion, scaled by the ratio of counts in the data to counts 
in the simulated spectrum. 

The proton yield is corrected for deadtime in the data 
acquisition system (10-20%, measured to better than 
0.1%) as well as several other small inefficiencies. Cor- 
rections for tracking efficiency, trigger efficiency, and par- 
ticle identification cuts (using the aerogel Cerenkov de- 
tector) were small (<2%) and independent of e. Finally, 
we required a single clean cluster of hits in each drift 
chamber plane to avoid events where the resolution is 
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FIG. 2: (Color online) Reduced cross sections as a function of 
e. The solid line is a linear fit to the reduced cross section, the 
dashed line shows the slope predicted from a global analysis of 
previous Rosenbluth results [jj , and the dotted line shows the 
slope predicted by the polarization transfer experiments 0. 



worsened by noise in the chambers. This significantly re- 
duced the non-gaussian tails in the reconstructed quan- 
tities, but led to an inefficiency of roughly 7%, with a 
small (0.25%) e dependence, possibly related to the small 
variation of rate with e. We corrected the yield for the 
observed inefficiency and applied a 100% uncertainty on 
the e-dependence of the correction. 

The absolute uncertainty on the extracted cross sec- 
tions is approximately 3%, dominated by the uncertainty 
in the angular acceptance (2%), the radiative corrections 
(1%), corrections for proton absorption in the target and 
detector stack (1%), the subtraction of endcap and inelas- 
tic processes (1%), and the uncertainty in the integrated 
luminosity (1%). We apply a tight cut on the solid angle, 
using only the data in the central 1.6 msr of the total «6 
msr acceptance. These tight cuts limit the elastic data 
to the region of 100% acceptance, but lead to the rela- 
tively large uncertainty in the size of the software-defined 
solid angle. Because the solid angle is identical for all e 
values at each Q 2 , this uncertainty affects the absolute 
cross section, but not the extraction of Ge /Gm ■ 

The largest point-to-point uncertainties, where the er- 
ror can differ at different e values, are related to the 
tracking efficiency (0.2%), uncertainty in the scattering 
angle (0.2%), the subtraction of the inelastic proton back- 
grounds (0.2%), and the radiative corrections (0.2%). 
The total point-to-point systematic uncertainty is 0.45%, 
and the typical statistical uncertainty varies from 0.25% 
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at Q 2 = 2.64 GeV 2 to 0.40% at Q 2 = 4.1 GcV 2 . The 
cross sections measured at 2.262 GeV have a slight addi- 
tional uncertainty (0.3%) because these data were taken 
at lower beam currents (30-50 /iA), and so are sensitive to 
any non-linearity in the beam current measurements and 
have somewhat different corrections for target heating. 

The reduced cross section, ctr = tG 2 m +eG%, is shown 
as a function of e in Fig. [2] The uncertainties shown 
are the statistical and point-to-point systematic uncer- 
tainties. Additional corrections, e.g. the effect of a fixed 
angular offset for all points, would lead to a change in 
the cross sections that would vary approximately lin- 
early with e. These uncertainties would not contribute 
to the scatter of the points or deviations from linearity in 
the reduced cross section, but can modify the extracted 
value of the slope. These uncertainties are dominated by 
uncertainties in the e dependence of the radiative cor- 
rections (0.3%), background subtraction, (0.25%), rate- 
dependent corrections in the tracking efficiency (0.25%), 
and the effect of a small beam energy or scattering angle 
offset (0.25%). These corrections yield a 0.55% uncer- 
tainty in the slope of the reduced cross section which is 
included in the uncertainties of the extracted form fac- 
tors. 

The form factors are obtained from a linear fit to the 
reduced cross sections. The results are given in Ta- 
ble [H] and the ratio h p Ge/Gm is shown in Fig. |3] to- 
gether with the results of previous Rosenbluth and po- 
larization transfer measurements. Note that for consis- 
tency with previous Rosenbluth measurements, the ef- 
fects of Coulomb distortion [l0( have not been included 
in the results in Table [D] Correcting for Coulomb distor- 
tion would lower H P G E /G M by 0.048, 0.042, and 0.032 
and increase G M /(p P G D ) by 0.009, 0.007, and 0.006 for 
Q 2 =2.64, 3.2, and 4.1 GeV 2 , respectively. 



TABLE II: Form factor values extracted from this measure- 
ment relative to the dipole form, Gd = 1/(1 + Q 2 /0.71) 2 . 



Q 2 


2.64 GeV 2 


3.20 GeV 2 


4.10 GeV 2 


Ge/Gd 
Gm/{^ p Gd) 
UpGe /Gm 


0.949 ± 0.040 
1.053 ±0.015 
0.902 ± 0.038 


1.007 ±0.052 
1.048 ±0.015 
0.961 ±0.051 


1.132 ±0.077 
1.031 ±0.015 
1.097 ± 0.077 



The results presented here are in good agreement with 
form factors extracted from previous cross section data, 
but have much smaller statistical and systematic uncer- 
tainties. This fact, combined with the consistency of the 
various Rosenbluth measurements , rules out most ex- 
planations of the discrepancy in terms of possible exper- 
imental error in the cross section measurements. The ac- 
curacy of the present work leaves little room for doubting 
that the Ge/Gm ratios reported from the Jefferson Lab 
polarization transfer experiments are inconsistent with 
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FIG. 3: (Color online) Extracted values of /j, p Ge/Gm from 
this work (circles), a global analysis of previous cross section 
data (Fig. 2 of Ref 01 1_ (crosses), and high-Q 2 polarization 
transfer measurements @,|6j (triangles). 



the form factors obtained using the Rosenbluth tech- 
nique, and makes it clear that the problem is not sim- 
ply experimental error in previous Rosenbluth measure- 
ments. The source of this discrepancy must be identified 
before the new insight into the proton structure provided 
by the recent polarization transfer data can be fully ac- 
cepted. 

One possible source for the difference between the two 
techniques is the effect of higher-order processes, or ra- 
diative corrections to the Born (one-photon exchange) 
cross section. The form factors are extracted from the 
cross section (or polarization transfer) measurements as- 
suming the Born approximation, so the effects of ad- 
ditional processes must be removed from the measured 
cross sections. We correct the data for higher-order elec- 
tromagnetic interactions such as bremsstrahlung, vertex 
corrections, and loop diagrams 8]. If we had detected 
the electron, the bremsstrahlung correction to the e de- 
pendence of <jr would have exceeded the e dependence 
coming from Ge- In this experiment, the e-dependent 
correction is much smaller and of the opposite sign. As 
the other radiative correction terms have almost no e 
dependence, the consistency between the new data and 
previous Rosenbluth results provides a significant verifi- 
cation of the validity of the standard radiative correction 
procedures. 

It has been suggested that higher-order processes such 
as two-photon exchange, not fully treated in standard ra- 
diative correction procedures, could explain the discrep- 
ancy 0, 0| . Such a correction would have to increase 
the cross section at large-e by roughly 6% relative to the 
low-e values. 0. This would mean that the Rosenbluth 
form factors, the form factors extracted using only the 
cross section data, would have large errors due to these 
missing corrections. While two-photon exchange will also 
affect the polarization transfer data, the corrections to 
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the polarization transfer form factors are expected to be 
smaller, though not necessarily negligible [IlLll3 |. 

Additional experimental and theoretical effort is nec- 
essary to determine if the discrepancy in the form factor 
measurements can be explained entirely by higher-order 
radiative corrections. The effect of multiple soft photon 
exchange (Coulomb distortion) has been examined |icj 
and yields a change in the slope of 1-2%, corresponding 
to a 3-5% reduction in the extracted value of /i p G e / 'G m , 
well below the level necessary to explain the discrepancy. 
Recent calculations [2 ^| show significant corrections 
to the cross sections due to two-photon exchange, though 
they appear to explain only half of the observed discrep- 
ancy. Additional effort is going into calculations of two- 
photon exchange for the cross section and polarization 
transfer measurements. It remains to be seen whether the 
effects of two-photon exchange are the reason for what 
must otherwise be considered a severe discrepancy. 

If missing radiative correction terms are responsible for 
the discrepancy, then the Rosenbluth form factors, and 
to a lesser extent the polarization transfer form factors, 
will not correspond to the true form factors of the proton. 
The discrepancy must be resolved before precise compar- 
isons can be made between models of proton structure 
and the measured form factors. While the Rosenbluth 
data may not provide the true form factors of the pro- 
ton, they still provide a useful parameterization of the 
electron-proton cross section, with the missing higher- 
order corrections absorbed into the extracted form fac- 
tors. Thus, the Rosenbluth form factors do reproduce 
the observed cross sections and provide the best parame- 
terization when elastic scattering is used to compare the 
normalization of different experiments, or when the elas- 
tic cross section is used as input to the analysis of exper- 
iments such as quasielastic A(e,e'p) scattering llflj. 

In conclusion, we have performed an improved Rosen- 
bluth extraction of the proton electromagnetic form fac- 
tors using detection of the struck proton rather than 
the scattered electron to decrease dramatically the un- 
certainty in the extraction. The results are as precise 
as recent polarization transfer measurements, but are in 



agreement with previous Rosenbluth separations and in- 
consistent with high-Q 2 polarization transfer data. The 
precision of these new results rules out experimental error 
in the Rosenbluth results as the source of the discrepancy 
between the two techniques, and provides a stringent test 
of the radiative corrections that are currently used in 
elastic e—p scattering. There are indications that this 
difference might come from two-photon exchange correc- 
tions, but we must better understand the discrepancy 
before precise knowledge of the proton form factors can 
be claimed. 
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